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ABSTRACT 



Context. Below 1 mHz, the power spectrum of helioseismic velocity measurements is dominated by the spectrum of convective 
motions (granulation and supergranulation) making it difficult to detect the low-order acoustic modes and the gravity modes. 
Aims. We want to better understand the behavior of solar granulation as a function of the observing height in the solar atmosphere 
and with magnetic activity during solar cycle 23. 

Methods. We analyze the Power Spectral Density (PSD) of eleven years of GOLF/SOHO velocity-time series using a Harvey-type 
model to characterize the properties of the convective motions in the solar oscillation power spectrum. We study then the evolution of 
the granulation with the altitude in the solar atmosphere and with the solar activity. 

Results. First, we show that the traditional use of a lorentzian profile to fit the envelope of the p modes is not well suitable for GOLF 
data. Indeed, to properly model the solar spectrum, we need a second lorentzian profile. Second, we show that the granulation clearly 
evolves with the height in the photosphere but does not present any significant variation with the activity cycle. 

Key words. Sun: activity - Sun: chromosphere - Sun: granulation - Sun: helioseismology - Sun: oscillations - Sun: photosphere 



1. Introduction 

The photosphere, the visible layer of the Sun, is the location 
where the energy transport previously dominated by convection 
and turbulence is largely done by radiation with an optical depth 
of 2/3. The gas is visible in the form of granules, that pene- 
trate inside the stable photosphere. These granules, as well as 
other larger structures like the mesogranules or the supergran- 
ules, are the manifestation of the different spatial scales of the 
convective motions occurring in this region o f the Sun (IZahnI 
119871; lRoudieretal]|1991l; lEspagnet etal]|1993h . 

The study of the granulation is particularly important in 
helioseismology because, on the one hand, it excites the so- 
called 5-min oscillations, i.e. the acoustic (p) modes, and, on 
the other hand, it dominates the power spectrum at low frequen- 
cies preventing the detection of low-order p modes. Indeed, in 
the case of velocity measurements, the lower de tection limit 
of ac oustic mode s is established ar ound 1 mHz dGarcfa et al.l 
1200 IL |2004a; Bro omhall et al. ] l2007h . To progress in the detec- 
tion of such modes as well as to increase the de t ection proba 



bility of gravity modes jAppourchaux et al 
2002HTurck-Chieze et al.l2004l;lGarcfa'erai 



20001: 
I2007t 



Gabriel et al 



Mathur et al 



2007h we need to better characterize the properties of the gran- 
ulation in order to reduce, if possible, their impact on the helio- 
seismic measurements. Th e coming GOLF-NG inst rument will 
soon address this problem dTurck-Chieze et ani2006l) . This new- 
generation instrument should improve the signal-to-noise (S/N) 
ratio of low-frequency modes by measuring the Doppler velocity 
at different heights in the solar atmosphere. It would thus benefit 
from the reduction i n the coherence of th e granulation with the 
atmospheric altitude dGarcfa et al ]l2004bh . 



In this paper, we analyze the Power Spectral Density (PSD) 
of velocity time sub-series from the GOLlQ instrument on board 
SOHC0 which is a solar disk integrated resonant spectrome- 
ter. With such an instrument, we can already study the mean 
behaviour of the solar granulation in some range of the at- 
mosphere. So t his wo rk comes in complement to the study of 
lEspagnet et al.l dl995h where the authors found that the photo- 
sphere is highly structured with two distinct layers below and 
above about 90 km . With GOLF, and the technique used in 
I Jimenez-Reyes et aP d2007h . we are able to study, without spa- 
tial resolution, a region located between 250 up to 550 km above 
the photosphere. We show that the granulation evolves with the 
height in the photosphere and that the granules tend to have 
shorter lifetimes with a weaker velocity when higher in the at- 
mosphere. 

Section |2] is devoted to the data analysis, with first a brief 
summary of the velocity calibration procedure of the GOLF sig- 
nal and secondly a description of the fitting procedure of the 
power spectra in more details. Section [3] is dedicated to a de- 
tailed study of the granulation motions and its evolution with 
time during the solar cycle. Finally, last section (Sect. [Hi will 
emphasize the main results of this paper and will anticipate on 
future incoming works. 



Global Oscillation at Low Frequencies jGabriel et aljl993) 
solar and Heliospheric Observatory jDomingo et al.|[T99^ 
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2. Data analysis 

2.1. Data 

GOLF is a resonant scattering spectrophotometer that measures 
the Doppler shift of the neutral sodium doublet (Di at A- 589.6 
nm and D2 at i = 589.0 nm). 



Height in the atmosphere 



inside 



V > 




Solar Na D 



Ib>Ih.o 

i,<i,„ 



outside 



GOLF referenced lines 



Fig. 1. Representation of the GOLF measurements. The Solar 
Na D line is displaced from the rest position due to a positive 
velocity (V > 0). The BW (blue wing) measures higher in the 
line which means deeper (close to the photosphere) in the solar 
atmosphere. The RW (red wing) measures closer to the bottom 
of the line, i.e. upper in the solar atmosphere. 



In the present study, we use time serie s of this experiment 
calibrated into velocity dGarcfa et al.ll2005T) . In this instrument, 
the light coming from the solar sodium absorption line (half- 
width 500 mA) traverses a sodium vapour cell, placed in a lon- 
gitudinal magnetic field of ~ 5000 Gauss, where it is absorbed 
and re-emitted in all directions. This scattered light is symmet- 
rically split into its Zeeman components allowing a measure- 
ment on either side on the wings of the solar absorption pro- 
file. The scattered photons are collected by two photomultiplier 
tubes. Due to a malfunction of the polarization subsystem that 
switches between the two wings of the sodium doublet, and to 
ensure a 100% of duty cycle, a single-wing working cycle was 
selected since April 11, 1996. Three different periods can be dis- 
tinguished during the operations of the instrument: from April 
1996 to June 25, 1998 GOLF has been observing on the blue 
wing of the sodium profile; after the SoHO recovery mission 
in September 1998, GOLF was restarted on the red wing up to 
November 18, 2002, when it was switched back to the blue-wing 
configuration. Since then, it has been running, unchanged, in this 
configuration. 

The main velocity field (+0.5 km/s) measured by GOLF is 
the SoHO orbital velocity. Due to this movement, the working 
points on the sodium profile change with a 1-year periodicity. 
Thus, GOLF is measuring at different altitudes in the solar at- 
mosphere (see Figure [TJ. Indeed when the orbital velocity Vorh 
increases, the measurements taken on the blue wing are closer to 
the photosphere, whereas on t he red wing the mea surements are 
further from the photosphere dGarcfa et al.ll2004bl) . 

The transformation from the observed counting rates to an 
intensity-like variable suitable for a velocity scaling includes two 
major steps: firstly, corrections to the raw measurements to re- 
move instrumental effects; secondly, compensation for the Sun- 
spacecraft variable distance. Then, before calibrating the data 
into velocity, corrections of some instrumental effects are ap- 
plied to the counting rates including those applied to the pho- 
tomultipliers, such as dead time, ageing of the phototubes and 
high voltage perturbations, as well as sodium cel l stem, pho- 
tocathodes, and filter temperature corrections (see iGarcfa et alj 



(l2005h for further details). In particular, in the second period of 
the blue-wing measurements, the correction of the photomulti- 
pliers temperature was no longer valid and it has been removed 
from the calibration pipeline. As a consequence, there is a small 
non-corrected offset in this last segment of data compared to the 
first period. Lastly, the calibration introduces a smooth high-pass 
filter with a cut-off frequency at about 1 x 10"^ Hz that produces 
an effect till the region of few days. To check the independence 
of our results from the calibration procedure, we have repeated 
the anal ysis using another velocity calibration procedure devel- 
oped bv luirich et all (|2000|) which involves a different set of cor- 
rections. The results were qualitatively the same using both cal- 
ibration methods. 

GOLF has been working for nearly 4200 days, a long period 
of time that allows a detailed study of the temporal evolution of 
its power spectrum during the whole solar cycle 23. To under- 
stand it and, in particular, the granulation, we are going to use 
two other quantities: 

- Th e MPSI (Magnetic Plage Strength Index ) , developed 
bvlChapman&BovdenI ( Il986h ; lUlrichI ( Il99ll) ; iParker et al.l 
(Il998h which can be used as a proxy for the mag- 
netic activity. This index is determined from the 150- 
foot solar tower magnetograms by measuring the sum 
(in all pixels) of the absolute value of the magnetic 
field strengths between 10 and 100 gauss. This number 
is then divided by the total number of pixels (regard- 
less of magnetic field strength) in the magnetogram (see 
http : //www . astro . ucla . edu/~obs/ 1 5S_data . html). 

- The GOLF observational height in the solar atmosphere 
computed from line-of-sight velocity taking into account the 
main displacement of the Na lin es from the rest position. We 
use the procedure described in iJimenez-Reves et al] (|2007|) 
to translate these velocities into a value of height in the solar 
atmosphere. 

It is important to remember that GOLF measures Doppler 
shifts of the Na line. Thus, it is mostly sensitive to a mixture of 
the vertical velocity towards the disk center and the horizontal 
velocity towards the solar limb. In fact, due to the single-wing 
configuration, the highest sensitivity will not be at the disk center 
but shifted towards the East or the West depending on the wing 
where the measurements are done (see Figure 3 of iGarcfa et al.l 

dMi). 

2.2. Fitting procedure 

The original GOLF time series are sampled every 20 s and we 
used series of 91.25 days shifted every 22.81 days. It is impor- 
tant to note that in the statistical analysis we will use indepen- 
dent subseries only, while in the figures we will plot all of them. 
This length of data has been chosen as a good trade-off between 
the frequency resolution and the number of independent series 
needed to correctly track the variations in the height of the atmo- 
sphere during a year Fast Fourier Transform (FFT) techniques 
were used to compute each individual spectrum. The normaliza- 
tion use d in the PSD is th e so-called "one-sided" power spectral 
density dPress et al.lll992l) . A first selection over the spectra was 
made. We considered only spectra with a duty cycle above 85%, 
and we eliminated 3 series in September, October and november 
2002 -just before coming back to the blue wing- because they 
contain a mix of both wings. 

The velocity power spectral density (PSD) can be explained 
by a model in which each source of solar convective motions 
is described by an empirical law initially proposed by Harvey 
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Fig. 2. Results of two different type of fits, as explained in the text, applied to a GOLF spectrum of an arbitrary taken subseries of 
91.25 days long. Top: Left, PSD with a fit using 8 parameters (one lorentzian) to adjust the p-mode envelope; Right, ratio between 
the PSD and the fit around the envelope of /9-modes. Bottom: Left, PSD with a fit using 1 1 parameters (two lorentzians) to adjust 
the /7-mode envelope; Right, ratio between the PSD and the fit around the envelope of p-modes. The dashed lines represent the 
limits inside which the fit is performed. The color used for the different fits are: gray for the super granulation contribution, magenta 
for the granulation contribution, blue for the p-mode envelope, green for the noise and red for the harvey function (the sum of the 
granulation and supergranulation contributions). 



(lHarvevlll98l iHarvev et al.|[T993l) and corresponding to an ex- 
ponentially decaying time function, plus one or two Lorentzian 
funct ions for the envelope of p-modes ( IVasquez Ramio et alj 
I2OO2I) and a constant for the photon noise. The fit of the spec- 
tra, including both non-periodic and periodic components, are 
thus expressed by: 



N A 2 

V 4o-fT; 



/■=i 



(1) 



where 



- P(v) is the power spectral density; 

- Nph is the photon noise in units of {m/s)^/Hz', 

- i corresponds to the non-periodic velocity fields; 

- i corresponds to the periodic components; 

- cr, and T, are respectively the rms-velocity {mj s) and the 
characteristic time {s) of the ;-th background component (the 
limit of the first sum varies depending on the number of 
non-periodic background components of the spectrum to fit); 



- Aj and vq^ are the power ((m/s) /Hz) and the central fre- 
quency (Hz) of the Lorentzian profiles to fit to the periodic 
components at the higher frequency region of the spectrum, 
while Fj sets its width (Hz). These M possible peaks to fit 
can be identified as the so called photospheric or/and the 
chromospheric component; 

- finally, cj (as well as bi) are decay rates. 



Before fitting the spectra, due to a large number of bins and 
because the points are compressed in a logarithmic scale at high 
frequencies, different local averages have been made until the 
data are equally spaced in the logarithmic scale. This allows to fit 
the model all along the frequency axis. The fitting is performed 
in logarithm of the PSD employing a st andard non-linear least- 
square method (Levensberg-Marquardt) (iPress et al.lll992l) . 

Two non periodic components are fitted for granulation and 
supergranulation with b, = 2. To fit correctly the p-mode enve- 
lope, it is necessary to use 2 Lorentzian profiles. In preliminary 
fits, Fj was obtained close to 1 x 10"^, then we fixed this co- 
efficient to this value for the rest of the fits. In previous work 
(Regulo et al. 2002), an additional factor (v/vq.)"' can be found 
in the Lorentzian term, but following these authors, we have de- 
cided to fix the coefficient a , to zero in our fits. 
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We have tried to isolate in the background spectrum different 
effects that are merged in the data: the dependence with solar ac- 
tivity and the dependence with the observing point in both wings 
(solar atmospheric depth). Which effect is dominant? In the fol- 
lowing, we note BWl the first period in the blue wing, RW the 
following red-wing configuration and BW2 the second period 
in the blue-wing mode. Thus, the fits have 1 1 coefficients, i.e. 
2 coefficients for the granulation as well as for the super gran- 
ulation, 3 coefficients for each of both Lorentzians used to fit 
the /7-mode envelope (the T factor being fixed), plus the photon 
noise. If we fit with 8 coefficients, only one Lorentzian is used 
for the acoustic-mode envelope. 

Figure |2] shows an example of spectrum where each com- 
ponent is represented by a different curve and the sum of these 
curves is the total fit (left-hand panels). The ratio between the 
PSD and the fit is also shown (right-hand panels). We see the 
difference between a fit with one Lorentzian (8 parameters, top 
panels) and a fit with two Lorentzians (11 parameters, bottom 
panels) for the adjustement of the p-mode envelope. These plots 
indicate the necessity of both Lorentzians to correctly adjust the 
p-mode envelope and thus decrease both excess power present 
near 0.003 and 0.005 Hz in the 8-parameter fit. In the following 
section, we will focus on the deep analysis of the behavior of 
the granulation noise. To do so, we model the PSD using only 
a fit with two Lorentzian profiles (the conclusions obtained for 
this convective pattern would be qualitatively the same when us- 
ing a fit with a single Lorentzian). The analysis of the other pa- 
rameters will be left for future incoming works. However, we 
note that the fit of the supergranulation is less robust than the fit 
of the granulation. The supergranulation region is certainly af- 
fected by the calibr ation of the data that filters this region slightly 
jGarcia et al.l2005b . and by the small number of points in this re- 
gion due to a small frequency resolution. If we used non-filtered 
series of 365 days, the fit of the supergranulation is more stable. 

3. Results 

3.1. Evolution of the granulation with time 

Figure |3] shows the temporal evolution of the fitting parame- 
ters for the granulation as well as the activity during the last 
solar cycle represented by the MPSI index. The rms velocity 
in each wing is different, with a lower value in the red wing. 
Besides, the velocity is higher in BW2 than BWl as a conse- 
quence of the change in the calibration of the second blue-wing 
period when the photocathode temperature correction were re- 
moved. Moreover we clearly see the signal dependence with or- 
bital velocity as a one-year modulation, enabling to probe dif- 
ferent heights in the atmosphere as we will see in the next sec- 
tion. This one-year modulation has an opposite phase between 
the blue and red wings, that can be clearly seen on the top pan- 
els. The blue wing is directly correlated with the orbital velocity 
(third pannel) while the red wing is anticorrelated. Moreover, the 
amplitude is larger in the red wing than in the blue one. 

3.2. Evolution of the granulation with the height and the solar 
activity 

The height in the photospher e is computed follow i ng the pro- 
cedure described in detail in iJimenez-Reves et aTl (l2007l) . and 
briefly summarized here: 

- The so - called response functions (hereafter RFs; see iMeinI 
([1971 : iBeckers & Mflkevl ([19751) ) of the Na Dl and D2 



Fraunhofer lines are computed: these functions measure the 
reaction of the line profile when the atmosphere is per- 
I xirbed locally at a g iven height. The quiet Sun model C 
dVernazza et al] [ 19811) was employed as a reference because 
it produces spectral line profiles (both photospheric and 
chromospheric) that are generally in good agreement with 
observations made at low spatial resolution. 

- For each line we computed the RFs to changes in velocity 
by applying smaU perturbations on the model at each optical 
depth. 

- Then, the RFs for the full line profiles were convolved with 
the wavelength filter response of the instrument. 

- The calculations took into account changes to the line-of- 
sight velocity at different times of year. 

Figure |4| shows the evolution of the granulation parameters 
the altitude in the solar atmosphere and the MPSI index. In or- 
der to provide a quantitative assessment of the observed changes 
for the velocity and the characteristic time, we performed a lin- 
ear regression with each of the parameters: the characteristic of 
each linear trend is summarized in Table[T]for the variation with 
height and in Table|2]for the variation with the MPSI index. We 
used data from independent series that did not overlap in time. 
This was done for each of t he blue- and red-wing period s follow- 
ing the same method than in lJimenez-Reves et al.l ( l2007h . For the 
evolution with the MPSI index, the trend has been calculated af- 
ter removing the trend with height. Errors issued from the gen- 
eral fit have been used to ponderate the linear fitting. For the 
clarity of the figures, only error bars for the independent series 
have been plotted. The best fit for each of the red and blue data 
sets are indicated by the black lines in the panels of Figure |4| 
We also computed the Pearson and Spearman-rank correlation 
coefficients for the data that are also presented in both tables. In 
Table[T] there are also coefficients for the linear regression with 
all wings: in the case of the velocity and in order to remove the 
offset, we have applied an adhoc correction simply by adding 
(respectively removing) the half of the difference in the velocity 
intercepts (i.e. (0.48 - 0.44)/2) to BWl (respectively BW2). 

The dependence with the depth in the solar atmosphere is 
well visible. As it was explained above, when the observations 
are done in the red wing, we are observing higher in the solar 
atmosphere than when the blue wing is used. We notice several 
interesting behaviours: 

- For each wing, there is a clear dependence with the height: a 
decrease of o-g,. and an increase of Tgr with the altitude in the 
photosphere, that can be interpreted as a longer lifetime for 
granules higher in the atmosphere. There is a strong correla- 
tion with height, more important for o-g,- than for Tg,- (fitted 
gradients significant at 4.4, 10.6, 4. 7cr and 2.8, 2.6, 1.8crfor 
each mode operation). For (Xgr, the difference in the intercept 
between each wing is significant, whereas the error bars for 
the intercepts for Tg,- are really too big to be significant. It is 
worth noting that the significance of the linear fit for r^r is 
more important when considering all confused wings. 

- Concerning the velocity, there is an offset between the two 
blue wings but the same slope and the same linear trend 
(slope and offset) for the characteristic time; this is normal 
because a difference in the temperature calibration between 
the two wings can affect cfgr but not the value of Tg,-. The 
offset in velocity is significant and certainly due to the dif- 
ference in calibration. 

- The difference of slope between the blue wing and the red 
wing for the velocity is weak. It seems that there is a slight 
break in the profil of o-gr with height, near 400 km. 
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Fig. 3. Temporal evolution of the granulation (cr and r, top and second panel respectively) of the GOLF PSD during the solar cycle 
represented by the MPSI index plotted in the bottom panel. The third panel represent the orbital velocity with time. Here, the vertical 
dotted lines indicate the date of the changes in the GOLF observing configuration. 

Table 1. Linear regression of the granulation with height. 



Epoch 


Parameter 


Intercept 


Gradient 


Significance 


x' 


















Of fitted gradient 












BWl 




125 ± 102 


0.91 ±0.32 


2.8cr 


1.66 


8 


0.71 


0.55 


1.60x10^' 




o-gr (mis) 


0.44 ± 0.02 


(-3.06±0.70)xl0-^ 


4.4(r 


1.99 


8 


0.79 


0.85 


6.53x10-^^ 


RW 


Tgr (S) 


238 ± 109 


0.59 ± 0.23 


2.60- 


0.94 


15 


0.62 


0.64 


1.03x10-- 




a-gr {mis) 


0.53 ± 0.02 


(-4.54±0.43)xl0-^ 


10.60- 


2.43 


15 


0.83 


0.79 


4.22x10-'' 


BW2 


Tgr (s) 


272 ± 82 


0.46 ± 0.26 


1.8cr 


0.81 


19 


0.39 


0.25 


3.04x10-' 




O-gr (mis) 


0.48 ± 0.02 


(-3.07±0.66)xl0-'* 


4.7o- 


1.50 


19 


0.70 


0.63 


4.12x10-' 


All wings 


Tgr {S) 


214 ± 24 


0.64 ± 0.07 


9.1cr 


0.93 


42 


0.84 


0.78 


1.49x10-^ 




O-gr (mis) 


0.47 ± 0.005 


(-3.34 + 0.13)xl0-^ 


25.7cr 


1.95 


42 


0.93 


0.88 


4.59x10-'^ 



° Actual numbers of independent data points employed in each of the regression analysis 
* Pearson correlation coefficient 

Spearman rank correlation coefficient 

Two-sided probability of Spearman coefficient 



- However, there is nothing significant with the MPSI index: 
the scatter of the residuals with the activity is too big. 



4. Conclusion and perspectives 

In this paper we have investigated the vertical structure and time 
evolution of the solar granulation by means of a novel methodol- 
ogy based on the analysis of the full-disk Sun-as-a-star Doppler 
velocity observations. Thus we have been able to study the ver- 



tical velocity fluctuations and lifetimes of the solar granulation. 
We have shown that the GOLF PSD can be correctly charac- 
terized by our model, a Harvey function with two Lorentzian 
profiles. 

This work extends the study of lEspagnet et al.l (1 1 9951) where 
they found that the photosphere is highly structured with two 
distinct layers below and above about 90 km. With GOLF, we 
study the photosphere above ^ 280 km and we showed that gran- 
ules tend to live longer with a weaker velocity when higher in 
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Fig. 4. Evolution of the granulation velocity and characteristic time versus the height (left) and the MPSI index (right). The linear 
trend with MPSI has been computed after removing the trend with height. + and □: respectively BWl and BW2; A: RW. The solid 
lines in each panel are the best fitting gradient of linear regressions for independent series and for the separate blue- and red-wing 
data sets. Note the jump between red and blue operating mode. All panels show overlapping data for illustrative purposes, with 
errors plotted only on the independant data. 



Table 2. Linear regression of the granulation (residuals after removing the trend with height) with MPSI index. 



Epoch 


Parameter 


Intercept 


Gradient 


Significance 


X~ 


















Of fitted gradient 












BWl 


r(5) 


-25 ± 17 


72 + 41 


1.8<T 


1.15 


8 


0.22 


0.19 


6.51x10-' 




tr {ml s) 


(-7.88 ± 4.09)xl0-^^ 


(2.25 ± 0.97)xl0-2 


2.3cr 


1.11 


8 


0.56 


0.60 


1.20x10-' 


RW 


r(s) 


-49 ± 43 


21 ± 18 


1.2o- 


0.83 


15 


0.04 


-0.02 


9.50x10-' 




tr (mis) 


(-19.14±7.39)xlO-3 


(0.82 ± 0.31)xlO-2 


2.60- 


1.88 


15 


0.28 


0.24 


3.98x10-' 


BW2 


T(S) 


-25 ± 15 


20 ± 15 


1.3<T 


0.58 


19 


0.45 


0.32 


1.88x10-' 




cr (mis) 


(7.96 ±3.76)x 10-3 


(-9.16±3.71)xlO-3 


2.5cT 


1.14 


19 


0.38 


0.41 


8.37x10-2 



° Actual numbers of independent data points employed in each of the regression analysis 
* Pearson correlation coefficient 

Spearman rank correlation coefficient 
'' Two-sided probability of Spearman coefficient 



the atmosphere. Following the results of lTitle et al.l (1 1 9891) . there 
is a strong correlation between the granule sizes and lifetime. 
Therefore, we can conclude that larger granules reach the top of 
the photosphere, while penetration of small granules decreases 
with the size. To be more precise, as Figure|4]shows a lifetime of 
about 400-550 s for granul es between 250 a nd 550 km, we can 
estimate from the article of lTitle et al.l (1 1 9891) and their figure 2 1 
that the granules in these altitudes have a lifetime-average size 
of about 1.2 arcsec and a maximum size of about 1.4-1.5 arcsec. 



We have found that the granulation rms velocities (cTgr) are 
between 0.28 and 0.4 m*-'. These values are very different from 
those obtained from high-resolution measurements which are 
about 1 kms^K The difference in spatial resolution is the most 
likely cause of this difference of three orders of magnitude. 

However we did not put in evidence a change with th e solar 
activit y. This is consistent with the result of limenez-Rev es et all 
( I2OO3I) that the solar cycle effects are very small compared to the 
change in the observing height in the photosphere due to the 
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orbital motion. It could be due to the fact that GOLF observes 
in the more stable part of the sodium lines. However if we do 
not find any signifi cant variation with the cycle, an other study 
jMuUer et al.l2007l) found a cyclic variation of the contrast of the 
granules, nearly in phase with the solar cycle, the contrast being 
smaller at the periods of solar maximum, but no corresponding 
variation in the scale. 

This work opens two perspectives: (1) to better understand 
the evolution of the whole solar convective background during 
the activity cycle, and specially the evolution of the acoustic- 
mode envelope presently characterized by the presence of two 
Lorentzians. We are likely to think that the excess power char- 
acterized by the second Lorentzian corresponds to the presence 
of chromospheric modes, (2) to extend this study to our new in- 
strument GOLF-NG observation. It will cover a larger part of the 
solar atmosphere thanks to the 8 extractions of the Doppler ve- 
locity with a proper determination of their location due to the 
measurement of both wings and the use of only the Dl line 
dJimenez-Reves et al.llIOOTflTurck-Chieze et al.ll2006i l2008h . 
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